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ABSTRACT 

At high angular frequencies, beyond the damping tail of the primary CMB power spectrum, the 
thermal Sunyaev-Zel'dovich (tSZ) effect constitutes the dominant signal in the CMB sky. The tSZ 
effect is caused by large scale pressure fluctuations in the baryonic distribution in the Universe 
so its statistical properties can provide estimates of corresponding properties of the projected 3D 
pressure fluctuations. The power spectrum of the tSZ is a sensitive probe of the amplitude of density 
fluctuations, and the bispectrum can be used to separate the bias associated with the pressure. The 
bispectrum is typically probed with its one-point real-space analogue, the skewness. In addition 
to the ordinary skewness the morphological properties, as probed by the well known Minkowski 
Functionals (MFs), also require the generalized one-point statistics, which at the lowest order are 
identical to the generalized skewness parameters. The concept of generalized skewness parameters 
can be further extended to define a set of three associated generalized skew-spectra. We use these 
skew-spectra to probe the morphology of the tSZ sky or the y-sky. We show how these power 
spectra can be recovered from the data in the presence of an arbitrary mask and noise templates 
using the well known Pseudo-C; (PCL) approach for arbitrary beam shape. We also employ an 
approach based on the halo model to compute the tSZ bispectrum. The bispectrum from each of 
these models is then used to construct the generalized skew-spectra. We consider the performance 
of an all-sky survey with Planck-type noise and compare the results against a noise-free ideal 
experiment using a range of smoothing angles. We find that the skew-spectra can be estimated 
with very high signal-to-noise ratio from future frequency cleaned tSZ maps that will be available 
from experiments such as Planck. This will allow their mode by mode estimation for a wide range 
of angular frequencies I and will help us to differentiate them from various other sources of non- 
Gaussianity. 

Key words: : Cosmology- Sunyaev Zeldovich Surveys - Methods; analytical, statistical, numerical 



1 INTRODUCTION 

Measurements of cosmological parameters from Cosmic Microwave Background (CMB) surveys such as those performed by satellite missions 
like WMAlQ and Planclfl can not only provide a very accurate picture of the background geometry and dynamics of the Universe but can also 
reveal much more detailed information about ongoing physical processes. Indeed, with all-sky coverage and wide range of previously uncharted 
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frequencies the data provided by experiments suc h as Planck will produce "secondaiy" science which is arguably as valuable as the "primary" 
science; see e.g. anim. Maiumdar. for a recent review. 

The large-scale properties of hot intergalactic gas can be probed through the multifrequency observations. Inverse-Com pton scattering of CMB 
photons kno w n as thermal Sunyaev-Zerdovich ef fect (tSZ) leaves a characteristic distortion pattern in the CMB spectrum lSunvaev & ZeldovichI 
h972lll98(]h : lBirkinshawl ( ll999t) : lRephaeiil ( Il995[) . The fluctuation of this distortion across the sky as probed by CMB observations can provide 
valuable clues to the fluctuations of the gas density and temperature. In the (low frequency) Rayleigh-Jeans (RJ) regime it produces constant 
decrement, while there is an increment at high frequencies; in between there is a a null (around 217GHz). This characteristic behavior is a 
potential tool for the separation of tSZ from the other contribitions to the temperature anisotropy. Based on accurate knowledge of the tSZ and CMB 
spectrum, foreground removal techniques have been developed to isolate the tSZ signal in the presence of primary anisotropy and instrumental 
noise. These techniques are extremely effective in the subtracti on of primary anisotropics due to its we l l understood (perfect black body) frequenc y 
dependence and almost exactly Gaussian statistical behavior l lLeachll2008l ; iBouchet & Gisper3l 19991 ; lOelabroullie, Cardoso & Patanchonll2003h . 
The tSZ effect is now routinely imaged in massive individual galaxy clusters where the temperature of the scattering medium can reach lOkeV. 
This in effect produces a change in CMB temperature of order ImK at RJ wavelengths. Individual galaxy cluster tSZ images have a variety of 
astrophysical and cosmological applications including direct measurement of the angular diameter distance to the cluster through a combined 
analysis of X-ray data and measurement of the gas mass which can be useful in estimation of baryon fraction of the Universe. The High Frequency 
Instrument (HFI) of Planck, in particular, has been designed with bands centered at the minimum, the null, and the maximum of the thermal SZ 
(tSZ) emission. The extraction of frequency cleaned CMB and tSZ maps of a catalog of galaxy clusters selected by their tSZ effect are part of 
the scientific program of Planck. Here we are however also interested in the general intergalactic medium (IGM) where the gas is expected to be 
at < IkeV in mild overdensites which leads to CMB contributions in the range. In this work we primarily focus on the statistical study of 
wide-field CMB data where tSZ effects lead to anisotropies in the temperature distribution due to both resolved and unresolved galaxy clusters, 
keeping in mind that the thermal tSZ contribution is the dominant signal beyond the damping tail of the primary anisotropy power spectrum. We 
primarily focus on analytical modelling morphological properties of the tSZ fluctuations using Minkowski Functionals (MFs). 

Modelling of the lower order s tatistics of the tSZ effect can be performed using either analytica l or numerical appro aches. Authors using 
analytical appr oaches (ISeliakll2000l;IZhang & Penll200ll; iKomatsu & Seii^l200ll ; IZhang & SettillioOTl ; ICooravll2000l 1200 1[) have gen erally used 
the halo model dCoorav & Sethll2002l)~ To trace the tSZ effect due to photo-ionized gas outside collapsed halos, see e.g. ICooravlJioOll) . A second- 
order perturbative formalism is used to model the bispectrum in this approach. The gas will typically be at a temperature similar to the ionization 
energy of hydrogen and helium. The bias associated with the pressure fluctuations is assumed to be redshift dependent. For tSZ effect from material 
within collapsed halos the shock-heated gas is typically assumed to be in hydrodynamic equilibrium in virialized h alos. The statistical descr iption 
of halos that include the number count distribution are ass umed to be described by Press -Schechter formalism dPress & Schechte3ll974) . The 
radial profile of such halos are assumed to be that of NFW jNavarro, Frenk & Whitelll996h . These ingredients are sufficient for modeling of tSZ 
effect from collapsed halos dCooravllioOol) . 

In addition to analytical modeling the numerical simulation of tSZ plays an important role in our understanding of the phys i cs involved 
jPersi et al.11199^; iRefregier et al.ll200d; ISeliak et all 1200 ll; ISpringe l et al 200 ll ; Iwhite, Hemguist & SDringelll200l iLin et al.ll2004l; IZhang et all 
l2004ICao, Liu & Fang|200HRoncarelli et alj2007l ; iHallman et al l2009l.i2007Flda Silva et all 19991) . Some of these studies incorporate complica- 
tions from additional gas physics such as radiative cooling, preheating and SN/AGN feedback, at least to a certain extent. The inputs are otherwise 
difficult to incorporate in any analytical calculations. On the other hand the simulations are limited in their dynamic range and can benefit from 
analytical insights. 

The tSZ power spectrum is known to be a sensitive probe of the amplitude of density fluctuations. Higher order statistics, such as the skewness 
or the bispectrum we study here, can provide independent estimates of the bias associated with the baryonic pressure. Typically a collapsed three- 
point statistics such as the (one-point) skewness is employed for this purpose. The skewnes s compresses all available i nformation in the bispectrum. 
The recently proposed skew-spectrum is a power spectrum associated with the bispectrum jMunshi & Heavensl2O10l) which is useful to probe non- 
Gaussianity as a function of scale or the harmonics I. In addition to the lower-order multispectra, morphological statistics suc h as the Minkowski 
Functionals (MFs) carry independent information of non-Ga ussianity and have been studied extensively in the literature dHikage et al.ll2008l 
l2006l ; iHikage. Taruva & Sutol '2003 ; Hi kage et al. I I2OO21 , 120031. At the lowest order the MFs are completely described by a set of three different 
skewness parameters that describe complete set of three Minkowski functionals in 2D. We extended these generalized skewness parameters to their 
corresponding skew-spectra. These skew-spectra sample the bispectra with varying weights and carry independent information. The estimation of 
these skew-spectra can be done at relatively modest computational cost. The skew-spectra can be constructed by cross-correlating suitable maps 
that are constructed from the original (frequency cleaned) tSZ maps. The construction involves differential operations on beam smoothed pixelised 
maps that can also be performed in the harmonic domain. The estimators for skew-spectra can be constructed following the general principle of 
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power spectrum estimation. We use tiie well known Pseudo C;s (PCL) approach developed by jHivon et alj|2002h . It can handle arbitrary sky 
coverage and arbitrary noise characteristics and beam patterns. We use the PCL approach to construct the variance in the estimators as well as to 
compute their cross-correlations. These results are accurate for surveys with all-sky coverage and can also be modified to take into account partial 
sky coverage, using the flat sky approximation. 

To relate to the e xperimental scenarios we will consider an experimental setup similar to the ongoing all-sky experiment Planck 
jThe Planck Collaboration 2006). We consider the range of harmonics (2, 2000) . However, the resu lts presented here are also applicable to 
smaller surveys such as Arcminute Cosmology Bolometer Array Receiver fACBAR. lRunvan et alj ( l2003h PI which covers the ^-range (2000, 3000). 
ACBAR is a multifrequency millimeter-wave receiver designed for observations of the cosmic microwave background and the Sunyaev-Zel'dovich 
effect. The ACBAR focal plane consists of a 16 pixel, background-limited, 240 mK bolometer array that can be configured to observe simultane- 
ously at 150, 220, 280, and 350 GHz with 4'-5' FWHM. Together with Planck these two experiments will cover the entire range of £ values up to 
3000. In addition to Planck we also consider a noise free ideal set up for the range of I values (2, 2000) as a reference. 

The paper is organized as follows. In ij2]we review the details of the analytical models involving the power spectrum and the bispectmm 
of the tSZ effect. In ij3]we present the skewness parameters that can be used to study the Minkowski Functionals. In ij4]we define the triplets of 
Generalized Minkowski Functionals. We propose the use of skew-spectra associated with individual generalised skewness parameters that carry 
more information than the ordinary MFs. We show how these skew-spectra are associated with the MFs and are related to the generalised skewness 
parameters. Next in !j5]we present the estimators that can be used to extract the skew-spectra from a noisy data set in the presence of observational 
mask and arbitrary beam. We also obtain the signal-to-noise for these skew-spectra for realistic observational scenarios. Finally, S|6]is devoted to a 
discussion of our results and future prospects. 

The particular cosmology that we will adopt for numerical study is specified by the following parameter values (to be introduced later): 
= 0.741, h = 0.72, fit = 0.044, D.cdm = 0.215, flu = + J^cdm, = 0.964, wo = -l,Wa = 0, as = 0.803, fl„ = 0. 



2 LOWER ORDER STATISTICS OF THE TSZ EFFECT 

In this section we briefly review th e two different approaches that are commonly used to model the tSZ effect (for more details see 
dCoorav. Baumann & SigurdsonlllOOSi) and the references there in). We will use these models later to study the morphological properties of the SZ 
effect. The tSZ temperature fluctuation 0^^(f2, u) = ST{Cl, u) /Tcmb is given by the (opacity-weighted) pressure fluctuations i.e. 

e^^(f^, v) = 5T{Cl)/TcMB = g.{x)y{Cl) = g,{x) / dr; a{f)) ar/m, 7r,(S^, j?); (1) 

Jo 

we will use the symbol x to represent a comoving coordinate. The electron pressure is denoted as tte = UekBTe. Here j/(n) is the Compton 
t/-parameter, ar is the Thomson cross-section, fes is the Boltzman constant, rUe is the electron rest mass, rie is the electron number density and 
Te the electron temperature. The function gv{x) encodes the frequency dependence of the tSZ anisotropics. It relates the temperature fluctuations 
at a frequency i/ with the Compton parameter y, i.e. g,^{x) — 2;coth(2;/2) — 4 and s = (/ii^/fcaTcMB) ~ i^/(56.84GHz). In the low frequency 
part of the spectrum gv{x) — —2, for x <ti 1, here x is the dimensionless frequency as defined above. Unless stated otherwise all our results 
will be for the low frequency limiting case. The conformal time rj can be expressed in terms of the cosmological density parameters Qm and 
by the following expression: ri{z) = dz' /}i{z') and E^(z) = H^(z)/Ho = [flmil + z)'^ + Qk{1 + z)'^ + Qa]- Introducing the r the 
Thompson optical depth and the fractional fluctuations in pressure as tt = Spe/pe we can write u) — Jg° drrTv{Cl, r). Here r is the 

comving distance, r is the Thomson optical depth, and overdots represent derivatives with respect to r. 

We will primarily be working in the Fourier domain. The three-dimensional electron pressure tt can be decomposed into its Fourier coefficients 
using the following convention that we will follow: 57r(x) = (7r(x) — (7r(x)))/{7r(x)) with (57r(k) = J d^x 57r(x) exp[— ik ■ x]. The projected 
statistics that we will consider can be related to the 3D statistics which are defined by the following expressions which specifies the power spectrum 
Pw{k), the bispectmm B-,r{k\, k2, fca) and the trispectrum r^(fci, k2, k^, ki) in terms of the Fourier coefficients: 

;ri,r2)(53D(ki +k2) (2) 
(,57r(ki,ri)57r(k2,r2)<57r(k3,r3))c = (27r)^B^(fei, fe, fca; ri)(53o(ki + ka + kg) (3) 
(57r(ki,ri)j7r(k2,r2)J7r(k3,r3)<57r(k4,r4))c = (27r)^T^(fci, fea, fe, A;4; ri)<53i3(ki + k2 + k3 + k4) (4) 

^ http://cosmology.berkeley.edu/group/swlh/acbar/ 



© 0000 RAS, MNRAS 000, 000-000 



4 Munshi et al. 



tSZ Power Sepctrum 



Planck 



+ 



2x10 



10" 





1000 2000 3000 

1 

Figure 1. The power spectrum of SZ is depicted for our choice of cosmological parameters. The halo model was used for the computation of the power spectrum. 
See section ^2.3l for a short discussion and relevant parameter values. For an all-sky Planck type experiment we consider the range / = (2, 2000) 



Our aim is to relate the statistics of pressure fluctuations (57r(x) to those of the underlying density fluctuations 5(x). The hierarchy of higher-order 
correlation functions can be defined using a similar set of equations. The corresponding multispectra will be denoted with a subscript 5 i.e. Bs 
and Ts will represent the hi- and trispectmm of the underlying density contrast and the power spectrum will be denoted as Pg. The S^d functions 
represent the 3D Dirac delta function and represent translational invariance of corresponding correlation hierarchy in real space. The unequal time 
correlators appearing in the expressions of the multispectra will typically collapse to equal time correlators due to the use of Limber approximations 
that we will discuss next. 



2.1 Linking projected statistics with 3D statistics 

The spherical harmonics decomposition of the projected field S^^(f2) will be represented as Of,n can be expressed in terms of a weig 

We will henceforth suppress the frequency v dependence of af^ as well as 6^^(f2). The line-of-sight integration over r projects the 3D (5;„i(k, r) 

onto the harmonics af^: 

efl = JdCl YiUCl) e^^(t7) = ][l f ° W^'^ir) J k^dkY,jiikr)YiUn)Si^{k,r); M/^^(r) = -26,(r)f; i = caTn,{z) (5) 

Here Yim{^) represents a spherical harmonic and ji {kr) is the spherical Bessel function of order /. The weight W^^ {r) is introduced by projecting 
the 3D statistics onto the sky. We use the following relation to relate (5;^ (k, r ) in the spherical basis with (5(k, r ) in Fourier basis to obtain that 

(5,™(k) = (27r)-'''' / dtl Yim{(l) 5(k) (6) 



. We can use this expression next in Eq.([5) and use the definition of power spectrum from Eq.(|2). Using the integral representation of the Dirac's 
delta function Ssd (k) = (27r) ~^ J exp(ik • x)d^x; in association with the Rayleigh's expansion of the plane wave in terms of spherical harmonics 
exp(ik • x) = 47r5]],„ i\ji(kx)Yim(0.k) Yimi^)- Usin g these expression we arrive at the following expression which projects the 3D power 
spectrum P{k,r\,r2) to the spherical sky ( ICooravir200lh : 



= - / W^^(ri)dri / W^"^(r2)dr2 / dkk'ji,{kn)ji^ikr2)PAk,ri,r2) (7) 



The following form of Limber's approximation dLimbeijl 954 ILo Verde & Afshordil20o"8h is remarkably accurate at large angular scales: I < 100: 
J F{k)k^dkji^{kir)ki2{k2r) — {Tv/2r^)F (l/r) Soiri — r2). The use of this form of Limber's approximation greatly simplifies the subsequent 
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Figure 2. The three one-point skewness parameters S^^^ (short-dashed), 5'^' (long-dashed) and S'^' (dot-dashed) are plotted as a functions of the FWHM 9 for a 
Gaussian beam. The results plotted are for ideal surveys without instrumental noise. The skewness parameter S^^"^ is the ordinary skewness describing the departure 
from non-Gaussianity at the lowest order The points are results from numeiical computations using the halo model and lines denote smooth fits to the results. The 
one point skewness parameters S''' is related to the skew-spectra 5^'*' through the following expression: 



results by reducing the correlator of multiple time slices to single time correlators dCooravlbOOlh . 







Here dA{r) is the angular diameter distance in terms of the comoving radial distance r. The Power spectrum Cf ^ is plotted in Fig-[T]as a function 
of the harmonics I. 

The bispectrum in the harmonic domain is represented by the following three-point correlation function: 



dSZ - f I I' ^" \ /oSZ^iSZ ^sz > 



(9) 



mm m 



The matrix above in Eq.([9} represents the Wigner 3j symbols jEdmondsIl 1 9681) and Bf^i^i^ defines the tSZ bispectrum in the spherical sky. The 
equation linking the projected bispectrum Bf^i^i^ and its 3D analogue B-n{k\, ^2, fca), defined i n Eq.([3j, can be derived following exactly the 
same procedure. We quote here the final expression which takes the following form ( ICooravir200lh : 



[W""" W] D f h h h \r _ i2h + l){2l2 + l){2h + l) 



Next we describe models for the 3D power spectrum P,r(fc; r) and bispectrum i3,r(fe, k' , k"; r) that will be used for computation of the projected 
SZ versions Cf^ and Bf^i^i^. Modelling of the tSZ statistics requires a model for the clustering of underlying dark matter distributions &. We will 
employ extensions of perturbative calculations as well as halo model predictions to model the underlying statistics of dark matter clustering. We 
will describe these model in following subsections. We have shown the angular power spectrum Cf ^ for SZ effect in Figure-[T]as a function of 
harmonic number I. For the computation of this power spectra a halo model approach was used which will be discussed in i]2.3l As we discussed 
in the Introduction, modelling the tSZ statistics involves modeling the underlying dark matter clustering and its relationship with the baryonic 
clustering and thence the fluctuation in baryonic pressure. This modeling has so far been performed using two complimentary approach. In the 
simpler approach the clustering of dark matter is described using second order perturbation theory or its extensions. The clustering of baryons are 
described using a biasing scheme. These inputs are sufficiently accurate especially at larger length scales to model the SZ sky. We will consider 
the linear biasing model first. 
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Figure 3. The absolute values of skew-spectra are plotted as a function of the harmonic I. From left to right the plots correspond to sf'"^ , S^"' and Sj^"^ respectively. 
We consider four different beams with FWHM = 0', 5', 15' and 25' respectively as indicated. These results correspond to ideal noise-free experiments and the 
resolution is taken to be Imax = 2000. The results are obtained using halo model prescription. 
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Figure 4. The absolute values of skew-spectra are plotted as a function of the harmonic I. From left to right the plots correspond to S'*^' , S'^' and 5^^' respectively. 
We consider four different beams with FWHM = 0', 5', 15' and 25' respectively as indicated. These results include Planck-type noise and correspond to an all-sky 
survey. 



2.2 Linear Biasing and Perturbative Approach 

In a perturbative approach the gravit ational dynamics is analyzed by expanding the large scale density field 5 in a perturbative series and a redshift 
dependent linear biasing is assumedtaoldberg & Spergel' 1999a'P). It is valid as long as the variance of the density contrast is smaller than unity 
and such a treatment is suitable for diffuse tSZ component (Hansen et al. 2005). For large smoothing scales such calculations can provide valuable 
insight to grav itational dynamics. The perturbative bispectrum Bs(k\, k2,kz; r) for the density contracst 5 from such calculations is given by 
jPeebIeslll97lh : 



Bs{ki,k2,kz\r) = F2(ki,k2)Ps{ki,r)Ps{k2,r) + cyc.perm. (11) 

2„-2/63 ki ■ k2 f kT_ k2\ 2 2/63 (k]_JC2)^ 

7 + 2fcifo {k2 ^ fcj + 7 " kfkl 



F2ikuk2) = i--n,/ +-^^1 + \+ n^/ ,.2.2 (12) 
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In the deeply nonlinear regime we do not have a complete analytical model of gravitational clustering. Typically, plausible approximations that 
include the halo model (to be discussed next) are used for this purpose, with varying degree of success. One such ansatz is the well known 
hierarchical ansatz tha t assumes that the higher order correlation functions can be constructed from the products of two point correlation functions 
teemardeau et all2002l) . The different tree diagrams thus constructed are all represented by various hierarchical amplitudes which are left arbitrary. 
The bispectrum in such a hierarchical scenario takes the following form: 

Bs{ki,k2, kr, r) = Q3 [P(fci)P(fc2) + cyc.perm.] ; Q-4n) = [4 - 2"]/[l + 2"+']. (13) 

The expression for Q-iiji) is adopted from ('Scoccimarro &Frieman"l999'). For a generic power spectrum one can replace n with the local linear 
power spectral index at (k\ + A:2 + A:3)/3 (Hui 1999). More elaborate schemes that interpolate the quasilinear regime and the highly non-linear 
regime have been devised in recent years Scoccimarro & Couchman (2001). Though we will restrict ourselves to the computation of the bispectrum 
in this paper, the approach can also be extended beyond third order. The hierarchical ansatz captures the salient features of non-linear clustering; 
the expressions are relatively simpler than the full perturbative expressions and have been tested in variety of cosmological contexts. In connecting 
the statistics of the electron pressure field tt namely the power spectrum P,r(fc; r) and B- ^lk, k' , k"; r) we w ill use the following local biasing 
model P^(k,r) = b^{rfPs(k,r) and B„{k,k',k") = b^(rf Bs{k,k,k;r). Following ^Goldberg & Sperge l 1999a.b; Cooray 2000, 2001) the 
bias is chosen to be of the following form h.„{r) = 6.^(0)/(l + z) with 6,r(0) = kBTe{0) / {meC^ bg) ■ These 3D Fourier statistics will be the 
ingredients for computation of projected statistics on the surface of the sky using 

EqM and Eq.flOll. 



2.3 Halo Model and tSZ: A brief Review 

The halo model description of large scale structure relies on modeling of the clustering of halos and predictions from perturbative calculation to 
model the non-linear correlation functions. In the context of halo model the dark matter is assumed to be in collapsed halos which are characterized 
by their mass, A brief Review radial profile, halo occupation numbers and underlying correlation hierarchy. The baryonic fluid is assumed to be in 
equilibrium with the da rk matter profile and i s characterized by a specific equation of state. The use of halo model is well established in cosmology 
and has a long history fcoorav & Sethll2002h . It has been used relatively recently is modeling of clustering of galaxies, weak lensing observables 
and statistics of CMB secondaries. In the context of tSZ it has been used to compute the leading order non-Gaussianity such as bispectrum as well 
as higher order statistics such as trispectrum which is important is computing the covariance of ordinary power spectrum. It is known that the halo 
overdensity at a give n position x, 5'^{x, M; z) can be related to the underlying density contrast (5(x, z) by a Taylor expansion as was shown by 
jMo&White1ll996h : 

5''(x,M;z) = [6i(A/;z)5(x,z)-|-i62(M,z)5'(x,2) + ...]. (14) 

The expansion coefficients are functions of the threshold v{M; z) — 5c{z) /a{M, z). Where Sc is the threshold for a spherical over-density to 
collapse and cr(M, z) is r.m.s fluctuation within a top-hat filter. The parameter 5c{z) is the value of the spherical over-density at which it collapses 
at a redshift z and is given by the following functional fit for a given Q,a and Q m '■ 



^ 3(12^ 
^ ' 20 



5 A , 1 - fij 



936 V O.M{l + zf 



(15) 



The rms fluctuation in the sphere of radius R is constructed from linearly extrapolated power spectrum P{k) g'^{M;z) = 
D^(z)/dlnfcA2(fc)|VK(fcP)p where W{x) is the to-phat window and A^(fc) is the dimensionless power spectrum is given by A^(fc) — 
k"^ P{k)/2-K^ . The linear order bias bi(M; z) introduced in in this expansion of 5h in terms of underlying density contrast S depends on the 
threshold iy{M; z) ( IMo. Jing & Whitelll997h : 

5c{z) d^{z){l + [av^{M-z)]i'} 

The linear growth factor D{z) is defined by the following expression: D{z) oc H{z) dz'{l + z')[H{z')]~^ and is normalized to unity at 2 = 
i.e -D(O) = 1. The main two ingredients in the halo model are the radial profile of the halos and the number density of halos. The number density 
is given by : f iu) = \/2A^a^ /nil + ( ai/^)~''] exp (— a!/^/2) and the threshold u. The parameters (p, a) takes the value (0, 1) for Press-Sechter 
mass function jPress & Schechtenl974l) . The constant A is determined by imposing the normalization J f{v)dv = 1. The halo model incorporates 
perturbative aspects of gravitational dynamics through the halo-halo correlation hierarchy. The nonlinear features take direct contribution from 
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the halo profi le; the number-density of haloes also needs to be determined. The total power spectrum -P*(fc) at non-linear scale can be written as 
( ISelia]dl200d) sum of two separate contributions to the power spectrum z) = Pl'^{k, z) + P^^{k, z). These represent contributions from 

clustering of dark matter particles in two different halos, i.e. the halo-halo term and a contribution from a single halo or one-halo term P^'\ These 
terms are in turn expressed in terms of two integrals: 

Pl^'iKz) = fiAKz); Pl\k,z) = [I^l'l{k,z)fP,{k,z) (17) 

The integrals themselves depend on the number counts of halos as a function of the concentration parameter c and the mass M. 



(18) 



4°](fc;^) = j dM j dc^£^JMk,\M,c;z)f; l[]l{k; z) = j dM j dc^^W{M, z)\-K,{k,\M,c- z)\ 

The bispectrum can be expressed in an analogus manner: 

Bl (fci , fca , fcs) = Bf (fci , fca , fcs) + (fci , fc2 , fca ) + B^'' (fci , fca , fcs) (19) 
Bl'' =l!{kuk2,ks); Bf (fci,fe2,fe3) =/2(fei,fe2)/?(fe3)B.(fc3) + cyc.perm. (20) 
Sf = [2J{ki,k2,hs)ll{k3) + lUk3)]Ii{ki)Iiik2)P6{ki)P^{k2) + cyc.perm. (21) 

The kernel J(fci, k2, fca) is same as F2{ki, k2, fca) introduced in Eq.(ll2t. We have introduced the following notation above: 

^ I'dM f dc^^b^iM,z)\7r4k,\M,c;z)\''. (22) 



The Fourier transform of the halo elctron pressure profile n{r) is denoted by ne{k) above: 



sin(fcr) 



TTe{k) = I d^r ne{r) exp(— ik ■ r) = / Anr^ dr ne{r) — ^ — -. (23) 







For the dark matter profile we assume pd{r) = PsU~^{i + with u = r/r^. Here the dark matter is assumed to be a NFW profile and 
is a characteristic scale radius. In the context of the spherical collapse model, the outer extent of a cluster is taken to be its virial radius 
r„ = [3M/4:npM{z)A{z)] where pm(z) is the mean background matter density of the Universe at redshift z and Am{z) is the over-density 
of the halo relative to the background density. The ratio of the virial radius to the scale radius is called the concentration parameter c = r„/rs. 
Together c and Af detemines the dark matter distribution of a given halo. The gas density profile in terms of the dimensionless parameter u is 
pg{u) = poPg{u). In this particular case a polytropic equation of state is typically used. According to the numerical simulation clusters of a 
given mass have a range of concentration parameters. The number density of clusters with a given mass M and concentration parameter c are 
expressed in terms of a bivariate distribution function d?n{M, c; z)/dMdc = dn/ dM{M\ z)P(c\M\ z). The conditional probability distribution 
function of a halo having concentration c for a given mass M in numerical simulations was found to be well approximated by a log-normal 
distribution. P{c\M\z)dc = exp [-(log c - log c)/2criog c] dc/(c Zn(lO)) with c{M,z) = co/1 + 2 [M/M*(z = 0)]""". Here Ah is the 
mass scale at which ^{M),; z) = 1 and Co = 9 and Uc = 0.13 and aiogc ~ 0.14. More complicated models are possible where a mass 
dependent aio„c can be incorpo r ated. F or more discussion on equilibrium baryonic gas profile and its relation with the underlying halo profile see 
ICoorav, Baumann & SigurdsonI ( l2005h . 



3 MINKOWSKI FUNCTIONALS 

In this section we will review the basics of the Minkowski Functional (MFs) in brief. The MFs are morphological descripto rs well known in 
cosmology. Put simply, morphological properties are the properties that remain invariant under rotation and translation; see e.g. lHadwigeJ 119591) 
for a more mathematical discussion. We will concentrate on a set of three MFs d efined in 2D that are defined over an excursion set E with a 
boundary 9S for a given threshold u. Following the notation of iHikage et alj j2008l) we define: 

Voiiy) = / da; Viiiy) = -[ dl; V2{v) f ICdl. (24) 

Je 4 Jg^ In Jg^ 

Here IC denotes the curvature along the boundary 9E. The three MFs Vo(/^), Vi{i') and ¥2(1^), correspond to the area of the excursion set 
S, the length of its boundary 9S as well as the integral of curvature IC along its boundary which is also related to the genus and hence the 
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Figure 5. The skew-spectra are plotted as a function of the harmonic I. From left to right the plots correspond to Si'^\ 5;^^' and 5p' respectively. We consider 
four different beams with FWHM = 0', 5', 15' and 25' respectively as indicated. These results correspond to an ideal all-sky noise fee experiment. A peiturbative 
calculation which assumes redshift dependent linear biasing was used. 
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Figure 6. The skew-spectra are plotted as a function of the harmonic I. From left to right the plots correspond to S'"', S^^^ and S^'^' respectively. We consider four 
different beams with FWHM = 0', 5', 15' and 25' respectively as indicated. These results include Planck-type noise. 



Euler characteristics. The MPs can be decomposed into two different contributions Gaussian (y') and non-Gaussian Wuiy) i.e. Vk(y) = 
^k'X) + 5Vk{y) with u — O/ao and (Jq — (O^)- We non-Gaussian contribution i.e. SVkif)- We will further separate out an ampli tude A in the 
expres sions of both of these contributions which depend only on the power spectrum of the perturbation through ao and ai (see e.g. lHikage et alJ 
( l2008l) ): 



^exp 



6 3 D 



A = 



Ul2 



(27r)(fe+i)/2 a;2-fca;fc \X2ao 



(25) 



(26) 



The constant cuk introduced above is the volume of the unit sphere in k dimensions. Wk = tt'°^^ /r(k/2 + 1) in 2D we will only need cjo = 1, 
uji — 2 and aj2 ~ tt. The Hermite polynoimals Hk{i') appear in the expression of MPs for weakly non-Gaussian random field. The lowest-order 
departure from Gaussianity 5vX\u) is encoded in three different generalized skewness parameters S'"', 5*'^', 5*'^'. The next to leading order 
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correction terms SVjf'^ (y) are typically neglected as they are order of magnitude smaller. Just as the generalized three skewness parameters they 
can be constructed from generalization of kurtosis. All three skewness parameters can be constructed from the bispectrum using different weights 
to sample different modes. In our next section we will discuss these generalized skewness parameters and introduce the power spectra associated 
with them. We will denote the power spectra with Sf'\ S^^"^ and S'P^ These will carry more information than their one-point counterparts the 
generalized skewness parameters. 

Estimation of MPs typically is done directly by employing a grid-based approach in real space, so that the integrals introduced in Ea.ll24t are 
estimated on a discrete set of points. Here we will show that the computation of 5*''' introduced in Eq.(|26|l can also be done in harmonic space. This 
will be particularly suitable for near all-sky surveys where galactic masks and point source masks can be incorporated in the analysis relatively 
easily. This will allow a cross-comparison of both methods for any systematics as well as residuals from component separation. We will deal with 
near all-sky surveys but flat-sky results can also be recovered using standard procedure which will be suitable for smaller patch-sky surveys. 

The Fig|2]shows the parameters 5*'"^ S^^^ and S^^^ as a function of FWHM assuming a Gaussian beam. We have assumed a halo model for 
computation of these generalized skewness parameters. The parameters defined in Eq. ( 126) . cro and ai, appear in the perturbative expansion of the 
MPs as well as in their normalization. These parameters depend on the level of noise. 



4 GENERALISED SKEW-SPECTRA AND MFS 

The skew-spectra are cubic statistics that are constructed by cross-correlating two different fields. One of the field used is a composite field typically 
a product of two maps either in its original form or constructed by means of relevant differential operations. The second field will typically be a 
single field but may be constructed by using differential operators on t he original fie ld. All of these three skewness contribute to the three MFs that 
we will consider in 2D. The first of the skew-spectra was studied by ( ICooravl l2001) and was later generalized bv Munshi & HeavensI (1201 01 and 
is related to commonly used skewness. The skewness in this case is constructed by cross-correlating the squared map [Q^ [Q)] with the original 
map [0(f2)]. The second skew-spectrum is constructed by cross-correlating the squared map against [V^0(f2)]. Analogously the third 

skew-spectrum represents the cross-spectra that can be constructed using [Ve{Q) ■ Ve(n)] and [V^e(f2)] maps. 

m ilia 

- .r ^ 2 y^|"^a+l)+/l(/l + l) + ^2(^2 + l)1giiii,Jiiii,M/ilVilTyi, (28) 
iDTTCTQCrJ •'y-' t J 

8^ [['(^ + 1) + + 1) - ^2(/2 + l)]Hl2 + 1) + cyc.perra^^Bu,lJu,l^WlWl,Wl^ 



(29) 



Jili2i3 _ /(2/l + l)(2/2 + 1) ( h h I 



13 



(30) 



^ ' " "2/3 + 1 V (2/3 + l)47r V 

=^(2Z + l)5'f' (31) 



i 



An 

I 



(32) 



This set of equations constitutes one of the main results of this paper The matrices here denote the Wigner 3j symbols lEdmondsl ( Il968h , Wi 
represents the smoothing window which can be, e.g., top-hat or Gaussian or compensated. Each of these spectra probes the same bispectrum 
Buii2 but with different weights for individual triplets of modes (I, hh) specifying a triangle in the haiTnonic domain. The skew-spectrum sums 
over all possible configuration of the bispectrum keeping one of its side I fixed. For each individual choice of I we can compute the skew-spectra 
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Sj^^ . The extraction of skew-spectra from data is relatively straight forward. It consists of construction of the relevant maps in real space either by 
algebraic or differential operation and then cross-correlating them in the multipole domain. Issues related to mask and noise will be dealt with in 
later sections. We will show that even in the presence of mask the computed skew spectra can be inverted to give a unbiased estimate of all-sky 
skew-spectra. Presence of noise will only affect the scatter. We have explicitly displayed the experimental beam bi in all our expressions. 

To derive the above expressions, we first express the spherical harmonic expansion of fields the [\7^Q{Cl)], [VO(f2) ■ VO(f2)] and [0^(f2)] in 
terms of the harmonics of the original fields Qim - These expressions involve the 3j functions as well as factors that depend on various k dependent 
weight factors. 



mi m2 — m 



[v"e(J7)];„ = J YiU^) [v'Qin)] = -i{i + 1)9;™ 

[e'(J^)]i,n = j dClYtUn) [e'(A)] = ^(-l)'"e;,„,e;,™, 

[ve(n) • ve(f2)];™ = f dnY;^(n)[ve(n) ■ ve(f2)] = ^ e,,„,e;,™, / dd y.uq) [vy,^,{{i) ■ vY;,™,(n)] (33) 

= IY1 + 1) + ^2(^2 + 1) - l{l + 1)] f dClYt*^{{l)Yi,,^, iCl)Yi,^, (Cl) 



h I2 I 



I Yi-^rMh + 1) + Hh + 1) - /(/ + i)]ei,^,ei„^2ihhi , 

; j m j 



(34) 



We can define the power spectrum associated with the MFs through the following third order expression: 

Vl'' = j:iV,Wl + 1) = i E(2^ + 1) \sl"'H{u) + U'^Hi.) + tit^^sl'^m^u) + ■■■}. (35) 



The "three-skewness" defines the triplets of Minkowski Functionals (MFs). At the level of two-point statistics, in the harmonic domain we have 
three power-spectra associated with Minkowski-Functional V^'"^^ that depend on the three skew-spectra we have defined. We will show later in this 
paper that the fourth order correction terms also have a similar form, but with an additional monopole contribution that can be computed from 
the lower order one-point terms such as the three skewness defined here. The result presented here is important and implies that we can study the 
contributions to each of the MFs Vk (v) as a function of harmonic mode I. This is especially significant result as various form of non-Gaussianity 
will have different I dependence and hence they can potentially be distinguished. The ordinary MFs add contributions from all individual I modes 
and hence have less power in differentiating various contributing sources of non-Gaussianity. This is one of main motivations for extending the 
concept of MFs (single numbers) to one-dimensional objects similar to power spectrum. It is worth pointing out that the skewness, along with the 
generalized skewness parameters, are typically less sensitive to the background cosmology. They are more sensitive to the underlying model of 
non-Gaussianity. The main dependence on cosmology typically results from the normalization coefficients such as ao and ai which are determined 
using the power spectrum of tSZ i.e. Q. 

In Figure-[3]we plot the three skew spectra for a range of Gaussian beam with FWHM = 0' — 25'. We consider an ideal all-sky experimental 
set up without any instrumental noise. The analytical results are computed using a halo model prediction for the bispectrum. The corresponding 
skew-spectra that include Planck type noise are displayed in Figure-]?] As expected the effect of noise is more pronounced at smaller angular scales. 
In Figure-]5]ideal noise free all-sky results are plotted for the diffuse component. The bispectrum in this case is constructed using the perturbative 
approach. The Figure-]6]shows the coiTesponding results for a Planck type all-sky experiments. Notice that for near all-sky experiments the signal- 
to-noise will degrade as fsky All results shown are for fsky = 1- 



5 ESTIMATORS AND THEIR SCATTER 

As we have noted before, the estimators for the skew-spectra can be most easily computed by cross-correlating maps in the harmonic domain. 
These maps are constructed in real space by applying various derivative operators. In the presence of mask the recovered skew-spectra will depend 
on the mask. The presence of mask typically introduces a mode-mode coupling. The approach we adopt here to reconstruct the unbiased power 
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spectra in the presence of mask exploits tiie so-called Pseudo-C; ( iHivon et alj|2002l) . which involves expressing the observed power spectra C; in 
the presence of a mask as a linear combination of unbiased all-sky power spectra. 

The three different generalized skew spectra that we have introduced here can be thought as cross-spectra of relevant fields. We denote these 
generic fields by A and B and will denote the generic skew-spectra as Sf'^. The skew-spectra that is recovered in the presence of mask will 
be represented as Sf'^ and the unbiased estimator will be denoted by Sf''^ . The skew-spectra recovered in the presence of mask Sf''^ will be 
biased. However to construct an unbiased estimat or Sf^'^ for the skew-spectra the following procedure is sufficient. The derivation follows the 



same arguments as detailed in iMunshi et all ( l2009h and will not be reproduced here. 



qA,B _ 

m I 



In this notation we can write the skewness parameters defined previously as: 

q(o)-Y qW-Y q(s>^^''©). c(2) - v ^(ve-ve.v^e) sy Y Y \-f— i 3_\ ni\ 

S =X(o)A, , b =X(,)S, , b =X(2)A, |X(o),X(i„X(2)| = |^^,^^^^,^j (37) 

The mode-mode coupling matrix i\f is constructed from the power spectra of the mask win and used for estimation of unbiased skew-spectra 
S^'^ . Typically the mask consists of bright stars and saturated spikes where no lensing measurements can be performed. The results that we 
present here are generic. The estimator thus constructed is an unbiased estimator. The computation of covariance of the the scatter in the estimates 
can be computed using analytical methods, thereby avoiding the need of expensive Monte-Carlo simulations. The scatter or covariance of the 
unbiased estimates {SSf'^ SSf!''^) is related to that of the direct estimates {SSf''^ Sbj^}'^) from the masked sky by a similarity transformation. 
The transformation is given by the same mode coupling matrix M. 

St'" = Y.^M-']u,st'''; (ssf'^'sst'^) = ^ A/,i^{5s^^^55f/^)M-,i,; (5^^) = s^^; {A, B} e {e,e^ (ve ■ ve), v^e} 

I' LL' 

(38) 

The power spectra associated with the MFs are linear combinations of the skew-spectra (see Ea.( l26t ). In our approach the power spectra associated 
with the MFs are secondary and can be constructed using the skew-spectra that are estimated directly from the data. 

The construction of an estimator is incomplete without an estimate of its variance. The variance or the scatter in certain situation is computed 
using Monte-Carlo (MC) simulations which are computationally expensive. In our approach, it is possible to compute the covariance of our 
estimates of various Sis, i.e. (SSiSSii) under the same simplifying assumptions that higher order correlation functions can be approximated as 
Gaussian. This allows us to express the eiTor covariance in terms of the relevant power spectra. The generic expression can be written as: 

IV1% = Et^^^'l"' (SVl'^SV^'^) = E MrL{S[V^%S[V^'\'}M£,\, (39) 

I' LL' 

We would like to point out here that, in case of limited sky coverage, it may not be possible to estimate the skew-spectra mode by mode as the 
mode coupling matrix may become singular and a broad binning of the spectra may be required. 



(40) 



Here the fraction of sky covered by the survey is denoted by faky The expressions for the skew-spectra are quoted in '^\s\"' 

and ®) given in Eq.( l32t . The expressions for covariance also depend on a set of power spectra i.e. '® ', S^'^ ^\ 

5,(ve ve,v 8) 5®'®. These are given by the following expression: 

Here Cf includes signal and noise power spectra Cf = Cfbf + . The experimental beam is denoted by hi. Using these equations it is possible to 
compute the scatter in various skew-spectra. The scatter will depend on the skew-spectra as well as various cross-spectra of various combination 
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Signal to Noise (S/N): Halo Model Results 




Figure 7. The signal-to-noise (S/N) squared for various skew-spectra S\ (left panel), Sj (middle panel) and S\ (right panel) are plotted as a function of angular 
harmonics I. We choose four different Gaussian beams with FWHM 6s = 0', 5', 15, 25'. The higher resolution (smaller FWHM) reaches higher signal-to-noise. 
The plots correspond to an ideal experiment with all-sky coverage and no instrumental noise. The expressions used for computation of the signal-to-noise are given 
in Ea. j45t -Eq. l47t . The results con'espond to an ideal all-sky fsky = 1 no-noise setup. For near all-sky esperiments the signal to noise will scale directly with fsky 



of product fields constructed from the original tSZ maps. 

_^ 1 r ro2 CA'i 



(42) 



^^^[ve.ve,v^e]^^[ve.ve,v^ei^ ^ /iy^ J^[v^e,v^e]^[v.v,v.v] ^ j^[ve.ve,v^e],2j ^44^ 
The cumulative signal to noise up to a given Z^ax using these expression for estimators S'"^ can now be expressed as: 

' = E (2' + l)(Sr^'^')V [x^oCP^'^^'cf'^' + [Sr'-'^'^f] (45) 

[(5/^)11] ' = f^^y E (2' + l)(5r''^''-")V [Xl,cf'^\l^'-^'^ + [Sf (46) 
;=2 

For the study of the tSZ effect we find that a robust determination of the generalised skew-spectra is possible for individual modes for almost the 
entire range of I values we have probed. This is important in differentiating them from various other sources of non-Gaussianity. The skew-spectra 
are integrated statistics, meaning that their value at a given / depends on the entire range of I values probed. These results can also be extended to 
take into account the cross-correlation among various skew-spectra extracted from the same data. 

^^5.[e^e]^^[ve.ve,v2e]^ = /-^ ^ ^gie^.v^B]^[B,ve-vB] _^ ^(e^ve-ve)^(e,v2e)j ^4^^ 

/rcie^.v^ejrcive-ve.v^e], _ 1 \^{e^,vB'VB)^(v^B,v^B) , ^[e^.v^e] ^[v^e.ve.vei] 

""^i / — /sky 21 I [ ' ' ~^ ' ' J ^ ' 

This discussion involves the lowest order departure from Gaussianity in MFs using a third order statistic, namely the bispectrum. The next to 



© 0000 RAS, MNRAS 000, 000-000 



14 Munshi et al. 



Signal to Noise (S/N): Halo Model Results 



1000 



3 100 



: S 



(0) 



\\\ 

I I 



1000 



100 



. (1) 



15' \ \- 

25' \ \e 

Mill I 1 1 1 1 1 1 1 1^ r 



II 1 1 1 1 1 

t- o(i) 



10 ^ s 
1 

0.1 
0.01 
0.001 



\1 

Vi 
I 



100 1000 

I 



100 1000 

I 



100 1000 

I 



Figure 8. Same as previous figure but with Planck noise added. Inclusion of noise degrades the signal-to-noise. The effect is most prominent for 5; which is 



(2) 

related to the fact that Sj gives more weights to smaller scales than other estimates and hence more affected by the noise. 

Signal to Noise (S/N): Perturbative Results 
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Figure 9. Same as Figure|8]but for perturbative calculations. The perturbative bispectrum has less power at higher I compared to halo model calculations. 
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leading descriptions are characterized by the trispectrum which is a fourth order statistics. It is possible to extend the definition of skew-spectra to 
the case of kurt-spectra or the power spectrum associated with tri-spectra. The power spectra associated with the Minkowski Functionals can be 
defined completely up to fourth order using the skew- and the kurt-spectra. However the corrections to leading order statistics from kurt-spectra 
are sub-dominant and leading order terms are enough to study the departure from Gaussinaity. In any case it is nevertheless straight forward to 
implement an estimator which which will estimate the power- spectrum associated with the MPs f rom a noisy data by including b oth third order 
and fourth order statistics is relatively straight forward. The issues has been dealt with in detail in ( iMunshi. Smidt & Cooravll2010l) in the context 
of CMB sky. 

In addition to the three generalized skew-spectra that define the MPs at lowest order in non-Gaussianity, it is possible also to construct 
additional skew-spectra that work with different set of weights. In principle an arbitrary number of such skew-spectra can be constructed though 
they will not have direct links with the morphological properties that we have focused on, in this paper they can still be used as a source of 
independent information on the bispectrum and can be used in principle to separate sources of non-Gaussianity, whether it be primordial or 
induced by late-time gravitational effects. In Pigure-|7]we have plotted the signal to noise for the three skew-spectra for different choice of the 
PWHM as indicated for an ideal experiment. The analytical results for the underlying bispectrum is obtained using halo model prescription. The 
left, middle and right panel correspond to S\'^\ S^' and S'P' respectively. In Pigure-[8]the signal-to-noise for an all-sky experiment with Planck 
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Signal to Noise (S/N): Perturbative Results 




100 1000 100 1000 100 1000 

I I I 

Figure 10. These results correspond to perturbative analysis but with instrumental noise for a Planck-type experiment included. 



type experiment is shown. As expected the with smaller beam size will have higher signal-to-noise. In Figure-|9]results for Planck type experiment 
is shown but using perturbative calculations and finally we show perturbative results that include noise in Figure-|9] The skew-spectra Sp' has 
higher signal-to-noise among the three skew-spectra. The very high signal to noise will allow an accurate determination of all the three spectra 
especially S;'"' and s\^^ for the entire range of / values probed. 



6 CONCLUSION 

In this study we have examined the prospects for extracting non-Gaussianity statistics from CMB surveys using MFs to probe the non-Gaussianity 
that originates from the tSZ effect. The tSZ effect is associated with the hot gas in large-scale structure that is detectable by multi-frequency 
experiments, in particular Planck.. When compared with the CMB temperature anisotropies the tSZ effect has a distinct spectral signature which 
means it can be effectively separated from the primary CMB contributions. This will provide an unique opportunity to probe tSZ effect using 
frequency-cleaned data in the very near future. In the past, the statistical analysis of the frequency cleaned tSZ maps has so far been mainly 
focused on lower order statistics. Morphological studies involving the MFs are complementary to the lower order statistics and carry independent 
information. In this paper we have generalized the study of MFs that are typically used in various other context in cosmology to the case of tSZ 
maps. 

As previously mentioned, the tSZ effect traces the fluctuations associated with the large scale distribution of the baryonic gas. The generation 
of pressure fluctuation in the virialized dark matter halos can be modeled by assuming them to be in hydrostatic equilibrium with the dark 
matter distribution within the halo. The shock heated gas typically correspond to over-densities in excess of 5 > 200. The unshocked photoionized 
baryonic gas typically traces the large scale distribution of the dark matter distribution. The typical over-density in unshocked photoionised baryons 
is typically 5 < 10. 

We have modelled the tSZ effect using two independent techniques. For smaller overdensities we employ a simplistic biasing model that rely 
on a perturbative description of dark matter clustering. The statistics of the hot gas is than linked to that of the dark matter using a linear (redshift 
dependent) biasing scheme. While this type of modelling is adequate for small overdensities a more elaborate analytical scheme is required for a 
detailed description of baryonic clustering at smaller scales. We consider a halo model based approach for the gas in collapsed virialised halos The 
specific number density and radial profile of these halos are modeled using Press-Sechter formalism or its variants. These two pictures of baryonic 
clustering is complementary to each other. The tSZ angular power spectrum corresponds to the projected power spectrum of the baryonic pressure 
fluctuation and the tSZ angular bispectrum correspond to the bispectrum of pressure fluctuations projected onto the surface of the sky. 

Few comments on the validity of the perturbative results are in order. The tSZ power spectrum by and large depends on the one halo term in 
the halo modelling. However in the perturbative regime we are mostly probing the Jean-Scale smoothed gas that is not in collapsed objects and 
traces the smoothed large scale dark matter distribution. To probe the large scale SZ effect removing X-ray bright clusters can reduce contributions 
from collapsed halos. Thus the effect captured by the linear biasing scheme should be understood as a signal in blank fields where such clusters are 
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absent. It provides the lower limit of SZ effect from large scale structure distribution. Tho ugh the diffuse com ponent of the tSZ effect is beyond 
WMAP detection threshold the situation may improve with future data sets such as Planck jHansen et al.ll2005h . 

The tSZ effect is intrinsically non-Gaussian. While the tSZ power spectrum is sensitive to the amplitude of the density fluctuations the higher 
order statistics such as skewness of tSZ effect can be used to separate the pressure bias from the amplitude of the density fluctuations. The skewness 
is related to the bispectrum of the tSZ effect. The individual modes of a specific bispectrum is defined by the triplets of the harmonic numbers 
{l-i , I2 , h). However individual m odes of the bispectrum has low signal to noise and may not be easy to estimate from a noisy data. In a recent 
study iMunshi & Heavens! ^20 id) advocated using a spectra correspond to the higher order spectra in general and bispectrum in particular. The 
skew spectrum is the lowest such spectrum that is constructed from bispectrum. The skew-spectra can retain some of the shape dependence in 
bispectrum without compressing it to a single number. Th e three different skew-spec tra that we introduced in this paper are generalizations of the 
ordinary skew-spectrum that was originally introduced in lMunshi & Heavens! ( |20I0 :1 and later used in the context of weak lensing and CMB for 
morphological analysis. These three different skew-spectra can also be used to study the morphology of the tSZ maps as described by the MFs. 
They associate varying weights to individual bispectrum modes when constructing the skew-spectra and carry independent information. 

Using near all-sky setup and noise that reflects ongoing CMB observations such as Planck we study how well the three different skew-spectra 
can be estimated from the data. We find that the data will allow a robust deteiTnination of two of the three skew-spectra that we have considered 
with a very high level of signal to noise. This is true for both perturbative results and the results that are based on halo model for the entire range 
of smoothing angular scales that we have studied. We also find the estimation of S'^' will be dominated by noise. The high signal-to-noise for 
the other two power spectra will allow mode by mode estimation of each skew-spectra. This can help to differentiate them from other sources of 
non-Gaussianity. 

The method that we pursue here depends on frequency cleaned tSZ maps. The tSZ effect can also be studied using cross-correlation techniques 
that involve external tracers. Such methods typically employ mixed bispectra. The results however lack the frequency infoiTnation and typically 
confusion noise dominated. The study of tSZ using bispectrum from frequncy cleaned maps typically provide additional signal-to-noise due to the 
frequency information. In the absence of frequency information the background CMB plays the role of intrinsic noise that degrades the signal-to- 
noise ratio. It is also interesting to note that removal of tSZ from the CMB maps may actually help in detection of other sub-dominant effects that 
we have not studied here, such as the kinetic Sunyaev Zeldovich effect (kSZ). 

Some of these techniques d escribed here will h ave wider applicability. The idea of generalized skew-spect ra was shown to be useful in th e 
context of weak lensing surveys jMunshi et alJ!201 ll) and probing primordial non-Gaussianity from CMB maps dMunshi. Smidt & Coorav!r201(]|) . 
In this study we have concentrated on only the leading order terms in the construction of the power spectra associated with the MFs. However 
the next to leading order terms can also be taken into account using the same formalism. The next to leading order terms will involve kurt-spectra 
that generalizes the concept of skew-spetra at fourth order. The four quadruplet of generalized kurt-spectra are related to the trispectra in a way 
that is similar to the relationship of skewspectra and the bispectrum we have considered in this paper. The generalized kurt-specta and the related 
generalized kurtosis can also be extracted from the data using the same PCL approach which we have discussed here. However the correction 
to generaised skew-spectra associated with MFs resulting from the trispctrum are expected to be negligible compared to the leading contribution 
from the generalized skew-spectra. 
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